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3. The jet in the first part of a swing is due to the local current created 
by the local difference of pressure; the vortex in the second part of the 
swing is due to the conjunction of the main stream with the opposing local 
current set up by the local pressure difference. 

My warm thanks are due to Mr, MacKinney for the zeal and ability he 
displayed in taking the instantaneous photographs from which figs. 3-6 are 
reproduced. 

DESCRIPTION OF PLATES. 

Fig. 3. — Obstacle with Pressure and Stream Indicators in Position in Tank, with the 

Water at rest. 
Fig. 4. — Swing from Left to Right—First Part. 
Fig. 5. — Swing from Left to Right — Second Part. 
Fig. 6. — Swing from Right to Left. 



Some Problems Illustrating the Forms of Nebulae. 

By George W. Walker, A.R.C.Sc, M.A., F.R.S., formerly Fellow of Trinity 

College, Cambridge. 

(Received March 27, 1915.) 

The possible forms of distribution of a mass of gaseous material under the 
influence of its own gravitation are of considerable interest in the nebular theory. 
The law of density which it appears most reasonable to assume is Boyle's 
Law, in which the pressure is proportional to the density, unless the pressure 
becomes so great that the material begins to resemble an incompressible 
substance. Although it is unlikely that the temperature is uniform through- 
out, still the solution under this restriction would be of value as a step in the 
direction of greater knowledge as regards possibilities in astronomical 
phenomena. 

The equations can be formed and lead to a differential equation for the 
surfaces of equal density. This equation is not linear, and in the three- 
dimensional case little progress to a general solution has been made. In the 
two-dimensional case, however, considerable progress can be made. 

A number of years ago I obtained the exact solution of the statical case 
of symmetry about the origin. Shortly after I found that Pockels had 
obtained the complete solution of the statical two-dimensional equation, and 
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Some Problems Illustrating the Forms of Nebulce. 411 

I drew attention to this in a short paper contributed to the 'Boltzmann 
Jubilee Volume/ p. 242, but, so far as I know, no further application of the 
solution obtained by Pockels has been made. I have recently examined some 
particular cases and have arrived at results indicating such a close analogy 
with actual astronomical forms that I venture to hope they may be of interest 
as suggesting types that may be expected in three dimensions. 

We begin by consideration of the fundamental equations. Let V be the 
gravitational potential, y the gravitational constant, p the pressure at any 
point of the gas, p the density at any point of the gas, x, y, the Cartesian 
co-ordinates of any point. 

The hydrostatic equations of equilibrium are 

1/dp dp\ /3V 3V\ 
p xdx' dyl \dx ' dy ) ' 

and the potential must further satisfy the equation 

a 2 v , d 2 v t 
w + df^ yp ' 

If p = pjh, where h is a constant, we get as an integral 

1/hAogp = V + constant, 
or p = poe hV . 

Hence Tf^+X 2 = —^irypo^. 

m 

This is the equation of which Pockels has obtained the solution in terms of 
two arbitrary functions as follows. As % = /i (# + iy) + /2 (# --- %) is the 

solution of -^ + tt^ = 0, then the solution for V and p is 

ox 2 oy J \ 

We require to put this in a form which gives real positive values to p, and 
I find that this can be secured by taking 

x=Mx+kj)+ M^¥Y 

Otherwise, if <£ and yfr are conjugate functions, so that 

we get y = + ^ + t J_., 

and hence p = P ^ = J^ {0^/3^ + 0^/3,^ } 



-iryh (02 + Tp + l) 



2 



412 



Mr. G. W. Walker. 



This remarkably general form gives us possible surfaces of equal density, 
when any known conjugate functions are assigned. I do not know whether 
a more general form giving real positive density exists, and it may be a matter 
of interest to pure mathematicians to investigate this point, just as it would 
be of value to applied mathematicians to know the most general form. 

Case I. — If we take 



(£> = (r J a) n cos nd , 



ty = {r/a) n sin n6, 



we get 



P = poe hY 



w 



a 



2 ' 



7ryha 2 ? 



2n 2 a? 



- — - sech 2 n losj (r/a) 



n 



iryha 2 7 



a ma 



,aj \r 



<a> n 



-2 



This is the most general case of circular symmetry. 

We note that when n < 1 the density becomes infinite at the origin, 
elsewhere finite, and vanishes at oo . When n — 1 the density is finite at 
the origin, and elsewhere is finite, vanishing at co . When n > 1 the density 
is zero at the origin, rises to a finite maximum, and then diminishes and 
vanishes at infinity. 

Writing the equation in the form 



p/pi = 4 



a' 



r 



n 



- + 



a 



,n 



-2 



the curves in fig. 1 show pj pi as a function of r/a, for the four special 
values n = 0, J, 1, and 2 ; a is still 'at our disposal while pi = 2n 2 jiryha 2 is 
determined when n and a are fixed. 
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Boyle's Law involves the possibility of an infinite density, and the law 
does not hold physically so far. In this and In other cases where a 
singularity occurs we may replace it by a solid nucleus arranged so as to 
give the proper value of V over its surface. 

This case of circular symmetry for n > 1 appears to have some analogy to 
a ring nebula. 

Case II — Take 

4> = h~ 1 \ogrja, i/r = b~ 1 ) 



2¥ 



,2 



we get p = P0 ^ = -^-j ± [log 2 (r/a) +ffi> + &»]-». 

The surfaces of equal density are given by 

P = -loo- 2 Z-W+ (J^Jf*. 

The curve, fig. 2, shows the particular surface 



2 — —log 2 1 + ~, 

a f 



drawn to a scale of 10. 




This case is somewhat peculiar. Starting with = when r/a = 1, the 
curve for positive values of proceeds as shown, continually approaching the 
origin by a succession of diminishing spirals. The locus is, however, 
symmetrical about = 0. If now we proceed to draw another surface of 
different density it would be found to cross the original curve at a succession 
of points for which > tt. This would mean that at such points two 
densities are possible, and we cannot admit this. We must, therefore, stop 
the curve at the nodal point = ir, and then the complete series of curves 
of equal density form a series of non-intersecting curves of this remarkable 
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pear shape, the density falling off to zero at infinity and increasing towards 
the origin with a singularity at the origin itself. The pear may be made 
sharper or blunter by taking different values of 6, which is at our disposal. 

This form of distribution is of some interest in connection with the- 
researches of Darwin, Poincare, and Jeans. 

Case III. — Let us take the elliptic co-ordinate transformation 

x-\-iy = c cosh (£4-^), 
and let <j> + iyfr = g-»(£-*o+*i). 

Then we get 

p = p e hY = - — — — (sin 2 77 + sinh 2 g)" 1 sech 2 ?i(£— £ )« 

2 Try he 

Thus the surfaces of equal density may be written in the form 

1 



sin 2 ?; = '° 1 



p cosh 2 n(t; — go) 
For n = 1, f = 0, we get 



l-cosh*f. 



sin 2 ?; = J ° 1 



1 



4-1 — cosh 2 £. 



p cosh 2 £ 

Fig. 3 shows the curves obtained for p/pi = -|, 1, 2. The density falls off 
to zero towards infinity, while for values p/pi > 1 the locus breaks up into* 

3 P/Pi 4 



FIG.3 



n= 1 




two oval curves about F and F', which are the foci of the original ellipses 
in the transformation. F and F 7 are singularities at which p would be 
infinite. The general character of the loci is not altered by giving f some 
finite value. For comparison fig. 4 shows on the same scale the result of 
putting n = 2. The equation is 

__ joi 1 



sin 2 ?; 



+ l — cosh 2 £, 



p cosh 2 2 £ 

and the curve drawn is for P / pi = 1. It is seen to be rather flatter than the 
corresponding curve for n = 1. 

This case then, as a whole, has some analogy in the case of a nebula with 
two fundamental equal nuclei. 
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Numerous other cases may be worked out, and it is simply a question of 
detail in assigning forms to the conjugate functions <f> and -v/r. But perhaps 
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these three cases are sufficient to indicate the generality of the method, and 
the remarkable interest of the forms that can be obtained. 

I pass to the consideration of an important related problem. 

The material of an actual nebula may be moving and not at rest, and in 
particular it may be rotating. As the rotation, if it exists at all, is very 
slow, the divergence from the statical equilibrium in such a case would be 
extremely small, so that the statical solutions are themselves of value. It 
is, however, of interest to know what the effect of motion may be. 

The simplest case we can consider is the final state in which the material 
rotates about an axis like a rigid body, with steady angular velocity a>. 
There the problem reduces to a statical one, when we take axes rotating 
with the material and add to the gravitational forces the centrifugal effect. 
In the case of a gas such a rotation, however small, can hardly be expected 
to extend indefinitely, and we should rather expect the motion to fall off 
and finally cease at a great distance. While keeping this in view, it is not 
without interest to examine what effect such an imposed uniform rotation 
would have. Mr. Jeans* has referred to this problem, but without giving 
details. Some of my results are in agreement with his conclusions and 
others apparently diverge. 

If the components of velocity are u and v, and we neglect viscosity, the 
equations for a steady motion referred to fixed axes are 

du . du 1 9 logo , 9V 

Ox oy ri ox ox 

dv . dv 1 9 log p . 9 V 

ox Oy h cy oy 

8,9 ex 

■^-pu + ^pv = 0. 
ex* oy 

9 2 V , 9 2 V 

* * Phil. Trans.,' vol. 213, p. 462 (1914). 
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In the case of symmetry, where p is a function of r only, we may take 

u = —yf> v = ^/, 

where /is a function of r only. 
The equations then reduce to 

p = p^F+ffrdr), 

1 3 3V 

-Trf-Tr = — 477*70. 

r dr or 

If / = a), we get the case of rigid body rotation, and p = / o ^ (V+ » a,2r2) . 
If / = cod 2 /r 2 , we get an ff irrotational " motion, and p = p ^^ v "i w3 " 4 /^). 

In the latter case the velocity falls off as r increases but is infinite at the 
origin. We might make / any function of r ; and among these we may 
select one, e.g. we~ r ^ a > which makes the velocity zero at the origin and also 
zero at great distances. But such cases are not steady and the influence of 
viscosity would arise. 

We cannot, in general, integrate in finite terms the equation 

r or or 

but if we suppose /so small that Y differs from its value V in the statical 
case / = by a small quantity x we have approximately 

O U'V 2 7h 

r ,r- r ^ = ~ — - — - . - (v + \f 2 r dr\ 

or or eosh J ?ilog(>/&) J 

while = A — =— — — : — — - — -— — {1 -f- A (v+ f/V dr) }. 

r 2 7T7Aa 2 r 2 cosh 2 ™ log (r/ a) l v * Jt/ y j 

Let %-f- J/V c?r = 2, and r = <ze*, then the equation may be transformed to 

3% _ 2w 2 2 3 20/>2 

« ■-- — — "~ r— — ~~Z ~r~ (X> ■~"i* 6 / . 

30 2 cosh 2 ?i0 30 y 

Now the equation 

3% 2n 2 A 

}- , ,% r=r (J 

3# 2 cosh 2 w# 

has a solution z = A tanh w#. 

Hence the equation may be written* 

coth nd 777: tanh 2 nd ^.z coth nO = a 2 ^ e 2d f 2 > 
60 06 06 

and the particular integral which we require is 

z == a 2 tanh ?&0 coth 2 «?i0 td,nhn6 ^f 2 e 2B d6 dd. 
* Cf. Boole, ' Differential Equations,' p. 205. 



Some Problems Illustrating the Forms of Nebulce. 417 

For uniform rotation/ = co, and then 

z = 2 a 2 m 2 tanh nd coth 2 n6 tanh n6e 2e d6 dQ. 



In the ease n = 1 this gives 



<X 2 &) 2 



3 = 



41og(y + l) + 



(y-i) 
(y+i) 



1 -,7 1 og(y+ 1 ) 



dy 



and for w 



<X 2 ft) 2 



^ = 



2 (tan 1 y- 



V 



1 + yV + $*+!) J (l-ytan- l yjrfy 



2 l \ 
where ?/ = r 2 /<x 2 . 

These curves are shown to scale in figs. 5 and 6, and we note that in both 
cases z = when r = 0,is everywhere positive, and finally becomes infinite and 
proportional to r^/a 2 when r is great. 



10 
9 



FIG.5 




The general effect of the rotation is then to increase the density from what 
would obtain in the statical case, and by increasing proportion as r increases 
from the origin, but even if it were permissible to pass to large values of r 
the density would still diminish to zero. We cannot, however, extend our 
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solution to large values of z, but we may say that the tendency of the rotation 
is (e.g. when n = 2) to make the density rise more steeply to its maximum 




value and to move the point of maximum density further out with increasing 
c» ; but there is no indication of a tendency to throw off* the outer layers. 
This latter result disagrees with the general conclusion stated by Jeans 
(loc. cit. ante). 

For " irrotational " motion / = o>a 2 fr 2 and the particular integral we 



require is 



z 



co 2 a 2 



.2y (y 2 +i) +tan y (y*+i)J - y 



y 



where y = r 2 ja % . 

In this case z is negative and infinite when r — 0, rises to a maximum 
positive value, and again declines to a negative and infinite value as r becomes 
indefinitely great. As before, we cannot extend our solution to great values 
of z, but we may say that the tendency of the motion is to make the distribu- 
tion hollow near the origin, to increase the concentration towards the 
maximum density at a finite value of r/a and finally to tend to throw off the 
outer layers when r becomes great. This latter conclusion is in agreement 
with Jeans. 

When we pass to other cases such as Cas^s II and III, where the original 
distribution of density is not symmetrical about the origin, it is rather 
difficult to think of any steady motion other than uniform rotation which 
would leave the distribution in anything like a permanent form. We must, 
therefore, confine our attention at present to this case, artificial though it 
may be. 
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Taking axes through the origin and rotating with the material with 
.angular velocity go, we have the equations 

p = ^(v+l-^ 

d 2 V d 2 V _ _ , 
dx 2 + dy 2 ~~ W * 

Again, let co be small and let V differ from its value V when co = by a 
;small quantity %. Then 

and ~^X + ~K = _ 4 7ry/oofe ^ v ° (% + 2 <° 2t2 )> 

ox* aif 

*or, if ^ + Jft) 2 r 2 = #, 

car dy J 

In Case II of the distribution of pear-shaped character I find that the 
.solution will have to be obtained by a somewhat complicated series, and as 
I have not yet gone very far with it I pass to Case III, which can be 
integrated in finite terms. We had as the statical solution 

p = p e hY ° = - — (sin 2 7] -j- sinh 2 f ) -1 sech 2 ntj. 

2 Try he 

Hence the equation we have to solve is 

(J Z f3 Z ^ 01 

^-ii + K~^ = — :=L Tr % (sin 2 ?? + sinh 2 ?) -1 sech 2 ^?-f 2co 2 , 

dx 2 oy 2 c 2 

or on changing to elliptic co-ordinates |, 77, we get 

5 2 + S= ^ T » + 6)V(cosh2f-cos2. 7 ). 

d£ 2 Crf COSh 2 w£ 

The solution is z = z 1 + z 2 cos 2r), where z x and z 2 are functions of £ only, 
and satisfy 

gg + _^i =w Vcosh2£ 
df J eostr ?££ 

■j 3% 2^ 2 (2cOSh 2 ^— 1) ■ 00 ra^onf^roKr 

..and -x— : - — - r-z — -z 2 = — w c respectively. 

a| 2 cosh 2 ?if F J 

A solution for z x , when 00 = 0, is z x = A tanh wf. Hence the particular 

integral is 

Zl = a> 2 c 2 tanh rcf J coth 2 «f J lanh raf cosh 2 £ d£ <Z£. 

A solution for z 2 , when co = 0, is £ 2 = A cosh 2 raf. Hence the particular 
integral is 



z 2 = ^~- H + i cosh 2 wf log cosh n%—±nt; tanh ?i| (2 cosh 2 ?i£ + 1) }. 

7h 2 
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The case n = 2 admits of easy expression, and we have 

ft) ' c 

Zi — —~ (cosh2£— 2 sech2£) 

ft) i 
z 2 = {| + -|-cosh 2 2f log cosh 2 £ — ^tanh 2f (2 cosh 2 2£+l)}. 

We may trace the effect of this on the density. Along the axis of x from 
to F (see fig. 4) £ = 0, and ?? changes from \ir to 0. 

Hence z = £1 + 22 cos 2 77 — — £o> 2 c 2 (l — -^cos 2 77), 

and so z ranges from ~-^ T co 2 e 2 at to -~^joo 2 e 2 at F. 

From F to 00 , ?; = and f ranges from to 00 . Hence 

2-= &1-4-22 = \to 2 (? {cosh 2^— 2sech 2 £ + -J- + -J- cosh 2 2 f log cosh 2 £ 

-i£tanh2f(2cosh 2 2£+l)} 

This function is — -/t a)2 ° 2 a ^ ^ increases and becomes positive as we increase 
£, and becomes infinite and positive with eo 2 e 2 /4cosh 2£ as f becomes infinite. 
Again, along the axis of y, r\ = \ ir, so that 
2 = z\ — z 2 = |&) 2 c 2 {cosh2| — 2sech2£~^— |-cosh 2 2£logcosh2£ 

+ i^tanh 2f(2 cosh 2 2f+l)}. 

At the origin £ = this function is — Y 7 T «¥, it increases as f increases 
and becomes positive, finally being infinite and positive with \ co 2 c 2 cosh 2 f as 
£ becomes infinite. 

The effect of the rotation is then to diminish the density near the origin 
and to increase it at greater distances, but there is no clear tendency to throw 
off the outer layers, since the function becomes positive and infinite and not 
negative and infinite. This is, again, in apparent conflict with Jeans' conclu- 
sion. It would thus appear from the cases considered here that a finite 
boundary of a gaseous mass is not to be looked for as a consequence of 
rigid body rotation, but rather as a consequence of some other type of 
motion in which viscosity may play a part. 

My sincere thanks are due to Sir Joseph Larmor for kindly criticism and 
advice. 

[Note added April, 1915.— It has been pointed out to me that in Case II 
dY/d6 becomes discontinuous along the line 6 = +7r. Hence this is a line 
of singularities. In accordance with the principle already stated we must, 
therefore, suppose this line to be occupied by solid matter of the proper 
density at any point, determined by the change of 1/r . dY/d0 in crossing 
this line. Theoretically this distribution would extend indefinitely, but 
practically the density of solid matter required at any point on 6 = it 
becomes very small as the distance from the origin increases.] 
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Fig. 3. — -Obstacle with Pressure and Stream Indicators in Position in Tank, with the 

"Water at rest. 
Fig. 4.— Swing from Left to Eight— First Part. 
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Fig. 5. — Swing from Left to Eight — Second Part. 
Fig. 6. — Swing from Right to Left. 



